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The influence of the size and shape of a dispersing and absorbing dielectric body on the local-field 
corrected spontaneous-decay of an excited atom embedded in the body is studied on the basis of 
the real-cavity model. By means of a Born expansion of the Green tensor of the system it is shown 
that to linear order in the susceptibility of the body the decay rate exactly follows Tomas's formula 
found for the special case of an atom at the center of a homogeneous dielectric sphere [Phys. Rev. 
A 63, 0538f 1 (2001)]. It is further shown that for an atom situated at the interior of an arbitrary 
dielectric body this formula remains valid beyond the linear order. The case of an atom embedded 
in a weakly polarizable sphere is discussed in detail. 

PACS numbers: 42.50.Ct, 42.50.Nn, 42.60.Da, 32.80.-t 



I. INTRODUCTION 

It has long been recognized that when an atom is 
situated in a material medium, the local electromag- 
netic field acting on it differs from the macroscopic one 
due to the gaps between the atom and the surrounding 
medium atoms, which are ignored on a coarse-grained 
macroscopic scale pi 0- Accounting for the difference 
between the two fields hence requires a correction — the 
local-field correction. Classical calculations of local-field 
effects, which typically have their origin in near dipole- 
dipole interactions, can be found in textbooks (see. e.g., 
Ref. |3|]). In quantum theory, investigations of local- field 
effects are often related to the problem of the sponta- 
neous decay of an excited guest atom embedded in a (di- 
electric) host. Local-field effects in spontaneous decay 
have been studied, e.g., for crystals and disordered 
dielectrics 0, 0, 0] on the basis on microscopic models 
for coupled atomic dipoles. 

In macroscopic descriptions, local-field effects are fre- 
quently taken into account by re gard ing the guest atom 
as being enclosed in a virtual 0, [l(| or real (spherical) 
cavity 0, El 0, El El E| surrounded by the medium, 
with the cavity size being small compared to the rel- 
evant transition wavelength. The cavity in the former 
model is virtual in the sense that it does not perturb the 
macroscopic field. It is filled by the atoms comprising 
the medium, which produce no net effect at the central 
position in important special cases such as cubic or ran- 
dom structures. In the latter model, the field is mod- 
ified by the presence of the cavity, which is an empty 
region containing only the guest atom. Microscopic mod- 
els often tend to agree with the virtual-cavity results 
0- HI* while many recent experiments on spon- 
taneous emission in dielectrics support the real-cavity 
model 0, El El El ■ It has been presumed that while 
the virtual-cavity model applies to interstitial atoms, the 



real-cavity model is specific to substitutional atoms, and 
that the case of substitutional atoms occursprevalently 
for impurity atoms in disordered dielectrics |5| ■ 

Since in the approaches to the local-field effects, the 
host medium has typically been assumed to be a bulk 
medium that extends homogeneously to infinity, the 
question of the effect of the size and shape of the host 
medium on the local field has arisen. In a macroscopic ap- 
proach, the spontaneous-decay rate of an excited atom in 
some free-space region can be given in terms of the imagi- 
nary part of the Green tensor of the macroscopic Maxwell 
equations, which characterizes the (macroscopic) envi- 
ronment of the atom. This relation in principle allows 
for including local-field corrections for atoms embedded 
in arbitrary material configurations by assuming a real 
(spherical) cavity surrounding the atom and calculating 
the corresponding Green tensor. Using the real-cavity 
model, Tomas [l3j studied the local-field correction to 
the spontaneous-decay rate of an excited atom, which is 
located at the center of a dispersing and absorbing dielec- 
tric sphere. Reformulating the result by representing it in 
a form, which does not explicitly refer to the highly sym- 
metric system considered, he made the conjecture that it 
may also remain valid beyond the specific example and 
hence also apply to other locations of the atom and other 
shapes of the host body. 

Based on a numerical computation of the respective 
Green tensors, Rahmani and Bryant 1141 ] considered the 
case of an atom at an arbitrary location within a dielec- 
tric sphere or a dielectric cube. Comparing their results 
for the dielectric sphere including the local-field correc- 
tion with earlier results disregarding the local-field cor- 
rection [53, EH- they suggested a rate formula, which 
in the case of weakly absorbing material corresponds to 
Tomas's formula. However, since their approach relies 
heavily on numerical calculations, it cannot produce ex- 
plicit expressions for the quantities that are related to 
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the local-field correction. 

The exact analytical evaluations of the Green tensors 
of realistic systems which have finite sizes and include 
a cavity can be very cumbersome. In this paper, we 
present an attempt to overcome this difficulty by writ- 
ing the Green tensor as a Born series in terms of the 
susceptibility, where in many situations one can restrict 
oneself to several leading-order terms. In particular, we 
show that to linear order the spontaneous-decay rate of 
a guest atom in a dielectric host body can be separated 
into a term representing the local-field correction to the 
decay rate in free space and a term related to the scat- 
tering Green tensor of the body without the atom — a 
result which exactly corresponds to Tomas's conjecture 
mentioned above. Furthermore, we show that for atoms 
which are situated at the interior of a macroscopic body 
Tomas's conjecture remains valid beyond the linear or- 
der. We illustrate the theory by discussing in detail the 
case of an atom embedded in a spherical dielectric body. 

The paper is organized as follows. The basic equations 
for the spontaneous-decay rate and the Born expansion 
of the Green tensor determining the rate are given in 
Sec. |n] They are used in Sec. IIIII to study the problem 
of the local-field corrected decay rate within the frame 
of the real-cavity model, and a proof of Tomas's formula 
is given. The examples of an atom embedded in a bulk 
dielectric medium or in a dielectric sphere are examined 
in Sec. II VI followed by a summary (Sec. IV|) . 



II. SPONTANEOUS-DECAY RATE 

Consider an excited two-level electric-dipole emitter, 
henceforth referred to as an atom, which is positioned at 
ta and surrounded by dispersing and absorbing dielectric 
bodies. The spontaneous-decay rate can be given in the 
form of HI E| 

2k 2 

r= -^d A -ImG(r A ,r A , WA )-dA, (1) 

he 

where d A and wa are the (real) dipole matrix element 
and (shifted) frequency of the relevant atomic transition, 
respectively, and &a = wa/c. The Green tensor of the 
bodies, G(r, r',w), satisfies the equation 



HG(r,r',uj) = 5(r-r')I, 

, ,2 



H = V x V x 



■e(r,w) 



(2) 
(3) 



(I, unit tensor) together with the boundary condition 

G(r,r',w)^0 for |r - r'| -»• oo, (4) 

where e(r, ui) is the frequency- and space-dependent com- 
plex permittivity which satisfies the Kramers-Kronig re- 
lations. Note that satisfaction of the boundary condit- 
ion 10} is ensured by assuming Im£r(r, uj) >0. 



Equation Q always applies when the atom is placed in 
some free-space region. Separating the Green tensor into 
bulk and scattering parts and respectively, 

G(r, r» = GW (r, r', u ) + G« (r, r', «), (5) 

and taking into account that in the case where the bulk 
part refers to free space, the relation 



I m G(°)(r A ,r A) u; A ) = ^J 



(6) 



holds (see, e.g., Ref. we may rewrite Eq. Q} as 

2k 2 

r = r + -^d A .ImG«(r A ,r A , WA ).d A , (7) 

fe 



where 



kid 2 



Tn = 



A "A 



(8) 



is the spontaneous-decay rate in free space. 

If the atom is embedded in a body, application of 
Eq. JJ} requires special care in two respects. Firstly, the 
coincidence limit of the bulk part of the Green tensor di- 
verges when the permittivity of the body is complex, as 
is the case in general. Only if material absorption can 
be neglected so that the permittivity can be regarded as 
being real, e(r A ,u; A ) ~Ree(r A ,a>A), this limit exists, 



ImG(°>(r A ,r A; 



to A) 



fe — I 

6ir 



(9) 



(see, e.g., Ref. |24(), leading to 
Ik 2 

T = x/ir + -^d A .ImG( 1 »(rA I r A , WA ).dA. (10) 

he Q 

Secondly, the Green tensor of the macroscopic Maxwell 
equations does not account for the fact that the local 
field felt by the atom is different from the macroscopic 
one in general. That is, even if absorption is neglected, 
the rate formula (|10|l is not complete because it lacks the 
local-field corrections. 



A. Born expansion 

To calculate the (scattering part of the) Green tensor 
for an arbitrary arrangement of dielectric bodies, it may 
be helpful to use an appropriate Born expansion. De- 
composing the permittivity as 

e(r,u) = e(r,w) +x(r,w), (11) 

and assuming that the solution G(r, r', u>) to the equation 



HG(r,r',u) = 6{r - r')I 



(12) 
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is known [where H is denned as in Eq. (0 with e instead 
of e] , the Green tensor can be written in the form of a 
Born series, 



[cf. Eq. (JBJ], and applying Eq. (|13|l . we may represent the 
scattering part of the Green tensor in the rate formula (Q 
as 



G(r, r', w) = G(r, r', + ^ A fe G(r, r', w) 

fc=i 

A fc G(r,r»=(-) J[ J tfsjxfa,") 
b'=i 

x G(r,si,w)-G(si,s 2 ,w)- • -G(s fe , r', w), 
as can be verified using the relationships 



^A 1 G(r,r',w) = ^x(r, w) [G(r, r', w ) 
-AxGCry.o;)], 
A fe G(r,r',u,) = ^(r,w)[A t -iG(r,r» 



(13) 



G (1) (r A ,r A ,w A ) = ^A fe G(r A ,r A , 



(20) 



fc=i 



(14) 



For r and r' belonging to a free-space region, substitution 
of Eq. ifTTjl into Eq. l(T4*| yields the first- and second-order 
terms in the Born expansion, AiG and A2G, respec- 
tively, as follows: 



AiG(r,r',w) = — - / d 3 sx(s, w) [oo'J - 06'u'u' 

lOTT J 



- a'buu + 6fe'(u-u')uu']e iVTT ^ {q+q ' ] 



(21) 



(15) 



[u = r s, g= fcti, a = a(q), b = b(q); u' = s — r', q' = ku' , 
a' = a(gO, 6' = &(</')], 



-A fe G(r,r»] for k > 1. (16) A 2 G(r,r» = _/ d 3 SX 2 (s,c)[ aa 'J- a&'u'u' 



The expansion (|13|l of the Green tensor is valid for ar- 
bitrarily spatially varying e(r, u>) and x( r j w )- Obviously, 
it may be very useful when \ can be regarded as being 
a (small) perturbation to e such that one makes only a 
small error by disregarding the higher-order terms. In 
particular, this is the case if the bodies are weakly polar- 
izable, as we shall assume in the following. 



B. Weakly polarizable bodies 

For weakly polarizable bodies, it is natural to regard 
the susceptibilities of the bodies as a small perturbation 
to the free-space permittivity e(r, 10) = 1 + irj (77 — *■ 0_|_), 
i.e., |%(r, u)\ <C 1. This means that we focus on frequen- 
cies that are sufficiently far from a resonance frequency 
of the dielectric material. Note that the small (positive) 
imaginary part of s ensures that G fulfills the boundary 
condition according to Eq. Q so that the spatial inte- 
grals in Eq. (fT4*|l converge. We have (see, e.g., Ref. Q) 



GUM-'-' l -^I+^{aI~buuy^\ (17) 



where 



r \ 1 i 1 > . / 1 3« 3 . _,. 
a = a(q) = - + _-_ b = b(q) = - + _-_, (18) 
q q z q a q q z q a 

u = r r', u = u/u, k — ut/c, q = ku. 

Separating the Green tensor into bulk and scattering 
parts in accordance with Eq. (JSJ), assuming the atom to 
be located in a free-space region such that 

toG(»'(r A ,r A , WA )=ImG(r A ,r A , WA ) = ^ I (19) 

D7T 



- a'buu + 66' (u ■ u')uu'] e lVT +^ {q+q ) 

+ el^/ d3six(si ' w) / d3s2x(s2 ' c " )eivTT ^ (<?1+912+92) 

x [aiai 2 a 2 I - aiai2b 2 u 2 U2 - ai6i 2 a 2 Ui2Ui2 

- biai 2 a 2 UiUi + ai6i 2 6 2 (ui 2 -u 2 )ui 2 u 2 

+ 6iai 2 6 2 (ui -u 2 )uiU2 + 6i6i2a 2 (ui -Ui2)uiUi2 

- &1&12&2(U1-U12)(U12-U 2 )U1U2] (22) 

[ui=r-si, Ui2=si-s 2 , u 2 = s 2 -r'; qi = kui, ai = a(qi), 
&,• = &(*) for i€ {1,2, 12}]. 

III. REAL-CAVITY MODEL 

Consider an excited two-level atom embedded in an ar- 
bitrary dispersing and absorbing dielectric body charac- 
terized by e(r, uj). In order to find the spontaneous-decay 
rate including local-field corrections, we employ the real- 
cavity model, that is to say, we assume that the atom 
is located at the center of an empty-space region of the 
form of a spherical cavity of radius Rq, as sketched in 
Fig. n^a) , where we have denoted the cavity volume by 
C, the volume of the host body without the cavity by B' 
(the overall volume of the host body being B = C U B'), 
and all the remaining space by V. Hence, the permittiv- 
ity of the system changes to 



£ioc(r, w) 



1 



if r € C, 



e(r,w) if r ^ C. 



(23) 



The cavity radius Rc is a model parameter representing 
an average distance from the atom to the nearest neigh- 
boring atoms constituting the host body; it has to be 
determined from other (preferably microscopic) calcula- 
tions or experiments. Note that the real-cavity model 
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B=C<JB' ^ 


V 





FIG. 1: Schematic illustration of the decomposition of (a) 
G*- 1 ' into (b) Gq and (c) Gg ; the dot indicates the position 
of the guest atom. 



is applicable provided that the unperturbed host body 
is homogeneous (and isotropic) in the region where the 
guest atom is implanted, 

e(r, lu) = e(r A , to) for r e C. (24) 



A. Linear approximation 

Restricting our attention to the first-order term in the 
Born expansion 1)20(1 . the (scattering) Green tensor cor- 
responding to £ioc(r,w) can be calculated from Eq. lj2T]l. 



where after some manipulations one obtains 

G« (r A , r A) w A ) = ^ / , *( s > ^A)e (2l ~" )9 

x [a 2 / + (6 2 - 2a6)uu] 

= Gg } (r A , r A , wa) + G<£ ] (r A , r A , w A ) (25) 

[u = r A — s, q = & A w], where, again to linear order in the 
susceptibility, 

Cx c (r A , r A , uix) - 167r2 

x / d 3 se (2i ^ )9 [a 2 J+(6 2 -2a&)uu] (26) 

JB'UV 

is the scattering Green tensor of a spherical cavity embed- 
ded in a bulk medium of susceptibility x( w a) = x( r A, w A ) 
[cf. Fig.QJb)], and 

G«(r A ,r A , WA ) = -^L^d 3 SX (s, WA ) e ( 2 ^ 

x [a 2 / + (6 2 - 2a6)uu] 

(27) 

is the scattering part of the Green tensor of the host 
body without the cavity [cf. Fig.^c)]. The result can be 
verified by applying Eq. J^J [together with Eqs. 1(13(1 and 
((21(1 ] to the two systems mentioned and using Eq. 1(241) . 

Equations lj2*5j) - l(2T|) show that in the rate for- 

mula (JZJ) can be written as the sum of two terms, where 
the first term, G^\ only depends on the cavity radius 
and the local permittivity of the host body at the po- 
sition of the atom, whereas the second term, is 
determined by the properties of the host body in the 
absence of the atom. Hence, the term in the decay rate 
which is proportional to Im Gq can be regarded as being 
the local-field correction to the uncorrected term propor- 
tional to Im G^ . The fact that the local-field correction 
additively enters the rate formula is due to the linear ex- 
pansion in X- Inspection of the second-order term in the 
Born expansion, Eq. ((22(1 . indicates that in general, terms 
depending on both the cavity and the (unperturbed) host 
body will appear which may lead to a breakdown of the 
additivity. However, as shown in Sec. IIII Bl there are 
situations where a generalization beyond the linear ap- 
proximation is possible. 

Equations l(26() and ((27(1 can be further evaluated by 
introducing a spherical coordinate system whose origin 
coincides with the position of the atom, 

d 3 s -f / d0 / sin 6» d6» / s 2 ds, (28) 

JO JO JRi(cf> : 6) 

where Ri(<j>, 0) and R (<fi, 9), respectively, refer to the in- 
ner and outer boundary areas of the integration volumes 
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sketched in Figs, ^b) and^c). Performing the radial 
integral in Eq. (|2(j|) . we find, on recalling Eq. i|18|) . 

G ( c\r Al r A ,uj A ) = lim lim F n (q) - F (q c ) = -F (q c ) 

(29) 



(qc = k A Rc) with 



F v [q(0,6)} = - 



16tt 2 



d0 / d6»sin(9<! e i2l - n)q 
(i Jo 



1 2i 5 i \ T / 1 2i 3 i , 
3?" 3?-3^ U+U -- + -^ |SS 



■4iEi[(2i-»7)«]( -/-si 



(30) 



[g((^>, 6) = k A R(cj),8); s — (cos ^ sin 0, sin sin 0, cos0); 
Ei(x), exponential integral]. Using the fact that qc is 
independent of and cf> as well as the relation 



'o Jo 
from Eq. I|30[) we obtain 



47T 

d0 sin 6* ss = — I, 



(31) 



F a (qc)=- k -^(^-^-^ +l )e^I. (32) 

In particular, when the radius of the cavity is much 
smaller than the atomic transition wavelength, 



q c = k A R c < 1, 
then Eq. (|32|l reduces to 



F (qc) 



&ax(^a) 
Gtt 



i 

9c 



7i 
6~ 



0(qc) 



(33) 



I. (34) 



Substitution of this result together with Eq. (|29|l into 
Eq. i|25|) reveals that to linear order in x, 



G< 1 >(r A ,r A ,WA) 



fcAX(^A) 
67T 
,(1) 



7?; 



(fc A i? c ) 3 &Rc 6 



G^(r A ,r A , WA ) + 0(fc A i? c ). (35) 



For a homogeneous host body which is star-shaped 
w.r.t. the position of the atom (i.e., every point on the 
outer boundary area can be connected to the atomic posi- 
tion by a straight line that lies entirely within the body), 

G^ as given by Eq. H27JI can be evaluated in a similar 
manner. Using again spherical coordinates and recalling 
Eq. i(TS|) , we may evaluate the radial integrals in Eq. l(2T|) 
to obtain 



G^(r A ,r A , w A ) = lim F v [q Q ( 
»j->o + 



-Fq(0) 



lim lim Fr,(q)-F o (0) = lim F v [q { 

r\— >0+ q— >oo J r\— »0+ 



(36) 



[q a = <Zo(0;#) = k A R o {(f),0)], where henceforth i? (</>, 0) 
denotes the outer boundaries of the host body. 

Substituting Eq. into Eq. , we find that to lin- 
ear order in x the local-field corrected spontaneous-decay 
rate of an atom within a body can be given as follows: 

r = r + r c + r< 3 1) , (37) 

where Tq is the free-space decay rate as given in Eq. |JSJ , 



re 

To 



Imx(^A) ImxK) 7Rex(^A) 
(fc A i? c )3 k A R c 6 



(38) 



and 



2fc| 
he 



d A • Im G$ (r A , r A , lu a ) • dp 



(39) 



Equations H37|) - (|39() are valid for an atom embedded in 
a weakly polarizable, dispersing and absorbing body of 
arbitrary size and shape. Note that only Tc takes into ac- 
count the effect of local-field correction, while which 
is simply determined by the scattering part of the Green 
tensor of the host body without the cavity, reflects the 
uncorrected influence of the size and shape of the body 
on the decay rate. In particular, the first two terms in 
Eq. H38fl obviously result from irreversible energy trans- 
fer from the atom to the surrounding matter. If these 
two terms are dominant over the last one, effectively no 
radiation is emitted. 

We conclude this section by making two remarks con- 
cerning the necessary conditions under which Eq. I|37|) 
provides a good approximation to the spontaneous-decay 
rate, (i) The permittivity often appears in the Green 
tensor as a common factor y/eujR /c [cf. Eq. I|82l) in 
Sec. HVB~T] . To be on the conservative side, one should 



then take x^aRd 



[with R n 



He, </>) R o(0,(t>)} 



rather than x as the small parameter, (ii) Equation l|38|) 
already contains terms of the order of Imx/(fcA-Rc) 3 - 
Thus Imx/(fcA-Rc) 3 , n °t Inix, should be much smaller 
than unity for the linear approximation to yield a good 
estimation of the decay rate. The cavity radius Rc 
represents the average distance between the atom and 
the constituents of the body. For, say, fc A i?c ^ 0.01, 
it is sufficient to require that Imx^ClO -6 , so that 
Imx/(fc A i?c) 3 «l. 



B. Beyond the linear approximation 

Tomas 0] has found that in the special case of an 
atom situated at the center of a homogenous, dielectric 
sphere, the real-cavity model (Fig. [21 in the case where 
l A = 0) leads to the relation 
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FIG. 2: Real-cavity model for an atom embedded in a homo- 
geneous dielectric sphere. 



G (1) ( 



r A ,r A ,w A J 



kA f3(e 



6tt L 2e + 1 (k A R c ) 3 



9(e-l)(4e+l) 1 



5(2e + l) 
3e 



k A R c 



+ i 



9e 5 / 2 
(2e + l) 2 



- 1 



2e + 1 



G^ 1) (r A ,r A ,w A )+0(A;A-Rc), (40) 



e= e(w A ), where Gg'^r',^ is the scattering Green 
tensor of the homogeneous dielectric sphere, whose geom- 
etry dependence is — for r = r' = r A at the sphere center — 
entirely given by its dependence on the sphere radius R. 
As already mentioned in Sec. he made the conjecture 
that this relation might be more generally valid for (i) 
bodies of arbitrary sizes and shapes and (ii) arbitrary 
positions of the atom, provided that the atom is not on 
the body surface. 

Let us study the validity of this conjecture in more 
detail. It can easily be seen that to linear order in x, 
Eq. (|40|l obviously reduces to Eq. I|35[) . so the results of 
Sec. IIII Al show that in this order the conjecture is true 
even under the more general conditions of inhomogeneous 
host bodies, provided that the requirement (Pljl for the 
applicability of the real-cavity model is valid. Moreover, 
we will demonstrate in the following that Eq. (|4U[) re- 
mains valid beyond the linear order in \, provided that 
the respective host body can be regarded as being homo- 
geneous in the vicinity of the atom, i.e., 

e(r,w A ) = e(r A ,w A ) for |r - r A | < (1 + v)Rc, (41) 

with v being some small positive number. 

We begin with the case of a homogeneous body, in 
which case satisfaction of the condition l|41l) simply en- 
sures that the entire cavity lies inside the body. We 
first recall that G^^(r A , r A , w A ) determines the electric 
field reaching the point r A , in which it has originated, 
after being scattered at the surfaces of inhomogeneity 
[cf. Fig. nja)]. At the cavity surface the electric field 
can be separated into two parts, namely, one part that is 



(multiply) reflected at this surface and eventually returns 
to the point r A , and one part that is eventually transmit- 
ted to the exterior of the cavity. The contribution due to 
the reflected part is obviou sly give n by Gq (r A ,r A ,w A ) 
[cf. Fig. ^b)] , which reads 



sly give] 

miii 



Gq (r A ,r A ,u A ) 



k A { 3(e-l) l_ 



6tt\ 2e+ 1 (k A Rcy 



9(£-l)(4e + l) 1 

5(2e + l) 2 
0(k A R c ). 



k A Rc 



+ i 



9e 5 / 2 
(2e + l) 



I 

(42) 



In an infinitely extended body [cf. Fig. ^b)], the trans- 
mitted part at a point r outside the cavity is determined 
by the Green tensor [2{| 



Gc(r, r A ,w A ) 



k A A 

47T 



(al — 6uu)e 



(43) 



[u = r r A , q = nk A u, n = *Je, and a and b according to 
Eq. (|18fl l. The coefficient A (which corresponds to A 1 ^ 
in Ref. [23) is given by 



A = 



h(z )[zoh^(zo)]'-[z j 1 (z )]'h^(z ) 



hiz^h^izi))' - £[z Ji(2 )]^i i; (zi) 



(44) 



where zq = k A Rc , z\ ~ nk A Rc (the primes denote deriva- 
tives w.r.t. zq and z\), and 



h{z) = 



sin z 



COM Z 



and 



(45) 



(46) 



are the spherical Bessel and Hankel functions, respec- 
tively Inspection of Eq. I|44|) shows that 



A = 



3e 



2e+ 1 



0(k A Rc). 



(47) 



Substitution of this result into Eq. (|43|) leads to 



3c 



G c (r,r A ,tj A ) = 2g "^ - G^ 0) (r,r A ,w A ) + 0{k A R c ), 

(48) 

where Gg ' (r, r A , lu a ) is the Green tensor of the infinite 
body without the cavity [according to Eqs. I|17|) and i|18|) 
with k 1— > nk]. In other words, the electric field trans- 
mitted through the cavity surface to a point outside the 
cavity is equal to the field that would be transmitted to 
the same point in the absence of the cavity, multiplied 
by a global factor. By means of the general symmetry 
property (see, e.g., Ref. |2J|) 



G(r,r» = G T (r',r,t 



(49) 
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we also have 



3e 



G c (r A ,r,w A ) = 2g "^ 1 G^ ) (r A ,r,uj A ) + 0(k A R c ), 

i.e., the electric field transmitted through the cavity sur- 
face from an arbitrary point outside the cavity is also 
equal to the corresponding result in the absence of the 
cavity, multiplied by the same factor. 

For a finite body, the electric field will be (multiply) 
reflected from the body's outer surface, eventually giving 
rise to a field at r A . Without the cavity, processes of 
this kind are taken into account by replacing the infinite- 
body Green tensor Gg ' with its finite-body counterpart 
Gb- Combining this observation with Eqs. (|48|l and (|50J) . 
we conclude, on recalling the linearity of Maxwell's equa- 
tions, that the electric field which is transmitted though 
the cavity surface, scattered at the outer body surface, 
and finally retransmitted into the cavity is given by 



3t 



2e + 1 



G<£\r A ,r A ,u; A ) + 0(k A R c ), 



(51) 



so combining Eqs. I|42l) and l|51[) . we arrive at Eq. Q40jl. 

In this derivation, we have disregarded processes in- 
volving scattering of the field at the cavity surface from 
the outside. Processes of this kind can indeed be ne- 
glected, because their contributions are of orders higher 
than 0(k A Rc). To see this, note that for a cavity in 
bulk material [Fig. IHb)] the field at a point r outside 
the cavity originating in a point r' outside the cavity via 
scattering at the cavity surface, is characterized by [2^| 



G«(r,r',a, A ) 



OO 

Y, [<M(r,r> A )+B£2V(ry, 



m— 1 



where 

(zo)[zijm(zi)}' - [Z03m{Z0)]'jm{zi) 



D 



N 



jm(zo)[2l^m (Zl)]' - [zoj m (z )]'h m ) (z 1 ) 
jm(z )[zijm(zi)}' ~ £[zoj m (z ))' jm(zi) 



j m (zo)[zih$ (zi)]' - e[z j m (z )]' 



(52) 

(53) 
(54) 



and the tensors M(r, r', oj a ) and ^V(r, r', uj a ) do not de- 
pend on i?o Using the relations 26] 



3m(z) : 
hgHz) 



(2m + 1) 
i(2m - 1)!! 



^[1 + ^z 2 )], 
l + 0(z 2 ) 



2^m-\-l 



(55) 
(56) 



one can easily show that 

B% = 0[(fc A i?c) 2m+3 ] , = 0[{k A R c f^ ^ 

so Eq. leads to 

G ( c\r,r',Lj A ) = 0[(k A Rcf]. (58) 



(57) 



Processes involving reflections of the field at the cavity 
surface from the outside hence start to contribute in third 
order of k A Rc and therefore do not need to be included 
in Eq. gOJ. 

So far we have demonstrated the validity of the rela- 
tion (|40|) for homogeneous dielectric bodies of arbitrary 
shapes provided that the atom is situated at the interior 
of the body, such that the condition l|41|) is satisfied. Note 
that this condition practically coincides with the condi- 
tion l|24|l for the applicability of the real-cavity model. 
The arguments given above can be extended to inhomo- 
geneous bodies. If the condition l|41l) is satisfied, one can 
divide such a body into a more or less small homoge- 
neous part containing the cavity plus an inhomogeneous 
rest. Equations iffify . ffiljl. and (|5%|) then still describe the 
propagation of the electric field inside the homogeneous 
part of the body, and the effect of the inhomogeneous 
part can be taken into account by the scattering at the 
(fictitious) surface dividing the two parts. Consequently, 
we are again left with Eq. 14UII . 

Substituting Eq. (gOl into Eq. and recalling Eq. JHJ), 
we can again cast the local-field corrected spontaneous- 
decay rate in the form of of Eq. (|37|l , where now 



To I 2e 



1) 



1 



9(e-l)(4e+l) 1 



1 (fc A i?c) 3 
9e 5 / 2 



5(2e+ l) 2 



and 



P (l) _ 



(2e + l) 2 



^id A -Im 

fen 



1 



k A Rc 
(59) 



3s- 



2e + l 

(60) 

[e = e(rA,<^A)]- Recall that Gg is the (unperturbed) 
scattering part of the Green tensor of the host body with- 
out the cavity. Needless to say that to linear order in x, 
Eqs. lO and (UjUJ) reduce to Eqs. (J2HI and (JSHJ), respec- 
tively. 

In particular in the case of weakly absorbing matter, 
it may be sufficient to retain in Eqs. (|59|l and (|60|l only 
terms to linear order in Ime. From Eqs. (|37|l . (|59|l . and 
(|60|l it then follows that when 



Imed A ■ ReG^\r A ,r A ,w A ) • d A 



Reed A • ImG B (rA,r A , u> A ) • d^ 



then r can be given in the form 



where 



3Ree 
2Ree + 1 



Ar, 



< 1 (61) 



(62) 



Ime 



Ar _ 

TV ~ (2 Re e + l) 2 (fc A i?c) 3 
9(14Rce + l) Ime 



5(2Ree + l) 3 k A R c 



(63) 
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and Tb is the uncorrected decay rate as given by Eq. i|10|) 

with GW(rA,rA,a;A) = (r^, ta, wa). Rahmani and 
Bryant |l4j concluded from a numerical computation of 
the Green tensor of a dielectric sphere that contains a 
(small) empty sphere at an arbitrary position inside the 
material that the local-field corrected spontaneous decay 
rate has the form T= / 2 Tb + ST, where the shift term 8T 
is due to absorption. Comparing with Eq. I|62[l [together 
with Eq. I|63[l ]. we see that this is indeed the case when 
the effect of absorption is sufficiently weak and in partic- 
ular, the inequality 1|61|) holds. However, from Eqs. (|37H . 
B59JI . and (jfi()|> it is clearly seen that in gene ral, T can- 
not be given in the form assumed in Rcf . |l4| . Note that 
already the analytical solution for the special case con- 
sidered in Ref. 13] implies that this form cannot be true 
in general. 



IV. EXAMPLES 



Atom in a bulk medium 



In the case of bulk material, we have G B X) = 0, so 
Eq. JSTJl simplifies to 



r = r n 



(64) 



where Tq/Tq is given by Eq. £3) GSG3, which to lin- 
ear order in \ reduces to Eq. (J3SJ. In particular, when 
material absorption can be neglected, e~Ree, we simply 
have [HI 



3c 



2e+ 1 



v^r 



(65) 



[which follows directly from Eq. 162fl with Tb = ^/eTo, 
cf. Eq. (| 1 (>|) ] . Note that in this case the virtual-cavity 
model leads to 



(66) 



(see, e.g., Ref. ^3)- ^ i s n °t difficult to see that to linear 
order in x both Eq. I|65|) and Eq. I|66|) lead to 



6 



r , 



(67) 



i.e., the first-order theory for the spontaneous decay rate 
does not distinguish between the virtual- and the real- 
cavity model |9|, provided that absorption can be disre- 
garded. 



B. Atom in a sphere 

Let us consider an atom in a homogeneous dielec- 
tric sphere as sketched in Fig. [3 and first calculate the 



spontaneous-decay rate to linear order in x- In particu- 
lar, we need to calculate as given by Eq. (JSHJ) with 
from Eq. (3EI together with Eq. (jSOJ). Noting that 



q (<j), 9) = k A l A cos<9 + ^K 2 - - cos 2 9) (68) 

(cf. Fig- EJ is independent of </>, performing the ^-integral, 
and using the fact that cIa ■ § = d A cos 9 for a radially (_L) 
oriented dipole and dA • § = <^a sin 9 cos <j> for a tangentially 
(||) oriented dipole, we find 



r (l)-L(ll) _ 
B 

where 



lim Im 

4 »;^o 



x(coa) / dxf r ^(q ,x 2 



(69) 



1 2i 5 i 
3q% 2,ql 3q a 2 
1 2i 3 i 
<lo 9o Qo 2 

u\^-z\Ei[(2i- V )q }, 



,(2i-T))q 



fl(q ,z)=f^[q ,(l-z)/2], 



(70) 



(71) 



and q according to Eq. JBSJl, with cos 9 being replaced 
by x. 

In order to exactly calculate [Eq. (|60() ]. we make 
use of the exact Green tensor for a dielectric sphere as 
given in Ref. (25|, leading to 



3r 0T [ 9ie 5 / 2 

°lml^ — £(2m+l) 



(2e + iy 



x m(m + 1)C, 



— - Im<^ — > (2m + 1 

v v ' m— 1 



A* 



(nk A l A ) 



1 2 



nk A l A 



(72) 



C . 3 , ( „,, (A) + C »(«^)'- 



(73) 



where 



ehWizJlzohWizo)]' - [z 1 h^(z 1 )]'h^\z ) 



(74) 



r 1 



j^OM^ (*>)]' - [*lj^(*l)]'ft#(*b) 



,(1) 



(!)/ 



(75) 



with zq = &A-R, 2i = nk A R. 
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1. Atom on center 

To compare the exact decay rate [Eq. I|37fl together 
with Eqs. ©, and (g>HJ) as well as Eqs. (J72J) and/or (j73jl ] 
with the decay rate obtained in the linear Born approxi- 
mation [Eq. together with Eqs. ©, PS)), and ], 
let us consider, for simplicity, the case where the atom is 
positioned at the center of the sphere. Setting Za = and 
hence, q Q — k A R in Eqs. l|^l) - (|7T|l . and carrying out the 
x integral in Eq. I|69|l . we derive 



r 



r 



(i) 



(76) 



with 



1 R 



giving 



-r lirn Im<^ x(w A ) 



2i 



k A R 



(fc A i?) 3 
i 

— I — 
2 



{k A RY 



(77) 



.(i) 



1 



2i 



(fc A i?) 3 (fc A i?) 2 



fc A i? 



,2ik A R 



(78) 



for any finite radius i?, and approaching zero in the limit 
R^oo (which has to be taken before performing the limit 
r] — * 0) . When material absorption is small enough such 
that Imx(w A )/(fc A i?c) 3 "C 7Rex(w A )/6, and the radius 
of the sphere is large k A R » 1, then Eq. (|3T|l together 
with Eqs. and (|75|) reduces to [x(w A ) — Re x( w a)] 



r = 



i 



7x(^a) 
6 



x(wa) 



cos(2fc A i?)r . (79) 



Note that Eq. l(79"l) differs from Eq. ljo7j) in the second 
term, which reflects the finite size of the sphere. 

When the atom is on center, only the terms m — 1 
in Eqs. and (JZHJ) contribute to the exact T^, 

cf. Eqs. I|55|) and (|56|l . and hence we find 



r Im 



9ie 5 / 2 



.(2e+l) 2 

Combining Eqs. (|37|) . (|59J) . and (|80|l . we arrive at 
'3(e-l) 1 9(e-l)(4e + l) 



(80) 



r = Im 



2e + l {k A R c ) 3 
9ze 5 / 2 



5(2e + l) 2 /c A i? c 



(2e + l) 



(i + cf) r„, 



which can be shown to agree with Eq. (|78|l to linear or- 
der in x- In a more sophisticated linearization of Eq. i|81|) 
with respect to \i it ma y be advantageous to leave the 
X-dependence in the exponentials appearing in the coeffi- 
cient Ci unchanged. In particular, for sufficiently small 



1.18 



u 



1.09 



li \\ '\\ / A > 


II 


/ V / \ 


I 


(a) 

1 



12 




1.15 - 



k A R 



FIG. 3: The local- field corrected spontaneous-decay rate 
(solid lines) and the corresponding result in the linear Born 
approximation (dashed lines) of an atom situated at the cen- 
ter of a dielectric sphere of permittivitiy (a) e = 1.1 + 10 -8 i 
and (b) e= 1.2 + 10 _8 i as a function of the sphere radius 7?. 
The radius of the real (empty-space) cavity is Rc =0.01/&a. 
For comparison, the decay rate in a bulk medium of the same 
permittivity (dotted lines) is also plotted. 



material absorption, Imx{u A ) / (k A Rc) 3 "C 7Re x(^a)/6, 
and for a sufficiently large sphere, k A R^> 1, this kind of 
linear approximation leads to [\(uj a ) — Rex( CJ A)] 



7x(^a) 
6 

x(^a) 



r„ 



cos [2Re(n)fc A i?]e" 



-2Im(n)fc A -RT 



(82) 



recall Eqs. and Comparison of Eq. with 

Eq. I|79|l indicates that with increasing size of the sphere, 
the validity of the linear Born approximation to the de- 
cay rate becomes less satisfying in two respects, (i) The 
oscillating term is not damped, and (ii) there is an ac- 
cumulated error in the phase. The former effect is in- 
significant as long as Im(x)fcA-R ^ 1- This restriction 
still allows for a large range of sphere sizes; for example, 
for Imx — 10 -8 and Aa = 2ir/k A ~ 1/im, the condition 
i?< 10 V m follows. 

Figure [3] shows the dependence on the sphere radius 



10 



1.5 - 




1.08 



10 7 ImX 



FIG. 4: The local-field corrected spontaneous-decay rate 
(solid lines) and the corresponding result in the linear Born 
approximation (dashed lines) of an atom situated at the cen- 
ter of a dielectric sphere of radius R = 2/&a as a function 
of Imx (Ree = 1.1), for real-cavity radii Rc = 0.01/Zca an d 
0.02/fc A . 



of the local-field corrected spontaneous-decay rates ac- 
cording to the exact equation (|81|) and the linear Born 
expansion [Eq. (|37|) together with Eqs. iJSSJl and (J7SJ] for 
two values of the permittivity. For the parameters used, 
we have Imx/ (^A-Rc) 3 = 10 -2 , so that the first two terms 
in the curly brackets in Eq. I|81|l , which arise from absorp- 
tion, are negligibly small compared to the last one. As 
expected from Eqs. i|79|) and i|82|) . the decay rate oscil- 
lates with the sphere radius around the bulk value. For 
small spheres, the linear Born expansion is seen to be in 
good agreement with the exact result. It is further seen 
that, as the sphere radius increases, an increasing phase 
shift develops between the exact rate and the rate in the 
linear Born approximation. A comparison of Figs. EJa) 
and (b) reveals that this phase shift is larger for larger 
values Rex- 

The dependence of the decay rate on the imaginary 
part of the susceptibility is illustrated in Fig. 0] for two 
values of the real-cavity radius. It can be seen that 
the disagreement between the curves representing the 
linear Born expansion and the exact result increases 
with increasing material absorption or decreasing cav- 
ity radius, i.e., with an increase in the combined fac- 
tor Imx/^A^c) 3 - In particular, if in the case where 
Rc = 0.01/fcA, Imx changes from 10~ 8 to 10 -6 , then 
Imx/(fcA-Rc) 3 changes from 10~ 2 to 1, and the agree- 
ment worsens from being very good to being moderately 
good. For the larger cavity radius Rc = 0.02/ k\, which 
physically means a more dilute medium, Imx/(fcA-Rc) 3 
varies from 10~ 3 to 10" 1 for the same variation of Imx- 
Throughout this range, the combined factor remains 
much smaller than one and a good agreement is observed. 




FIG. 5: Position dependence of the local-field corrected (up- 
per curves) and uncorrected (lower curves) spontaneous-decay 
rate in the linear Born approximation of an atom in a di- 
electric sphere of permittivity e = 1.1 + 10 _8 i and radius 
(a) R=l/k\ and (b) R = 5/k\- The real-cavity radius is 
Rc = 0.01/fcA- The solid (dashed) curves refer to radially 
(tangentially) oriented transition dipole moments. 



2. Atom off center 

Basing on a numerical computation of rj 3 1 - ) " L( -"' ) as given 
by Eq. 169|l in the linear Born approximation, we have 
also studied the case when the atom is localized at an 
arbitrary position inside the sphere. Figure [3] compares 
the position dependences of the local-field corrected de- 
cay rate, as given by Eq. l|57|) together with Eqs. 
and l|fj§)) , and the uncorrected rate according to Eq. l(TT)|) 
(which is valid if absorption is neglected), with the second 
term being given by Eq. I|69f) . for two sphere radii (for the 
uncorrected rate beyond the linear Born approximation, 
see also Ref. H3]). From the figure it is seen that when 
I A > 0, i.e., when the atom is not localized at the center 
of the sphere, radial and tangential dipole orientations 
have to be distinguished, especially when the sphere ra- 
dius exceeds the wavelength of the radiation emitted by 
the atom, such that interference effects begin to play a 
role. Note that the whispering gallery resonances which 
may give rise to strong enhancement near the rim of the 
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FIG. 6: Decay rate of an atom located at (a) I a = 1/&a an d 
(b) /a = 5//A;a off center of a dielectric sphere of permittivity 
£ = 1.1 + 10~ 8 i as a function of the sphere radius in linear 
Born approximation. The cavity radius is IcaRc =0.01. The 
solid (dashed) curves refer to radially (tangentially) oriented 
transition dipole moments. For comparison, the curves for an 
atom at the center &a'a = are also shown. 



sphere are not manifested here. The existence of these 
resonances requires larger values of the real part of the 
permittivity or larger sphere radii. Unfortunately, in such 
cases the linear Born expansion is too rough an approx- 
imation to the decay rate. From the figure it is further 
seen that for the parameters used, the local-field correc- 
tion increases the decay rate; the amount of increase is 
determined essentially by the last term in Eq. 

In Fig. the dependence of the local-field corrected 
decay rate on the sphere radius is illustrated. It can 
be seen that the further the atom is displaced from the 
sphere center, the smaller the amplitudes of oscillation of 
the decay rate become. 



macroscopic Maxwell equations and expanding the Green 
tensor in a Born series enables one to systematically ap- 
proach arbitrary geometries. In particular, in the case of 
weakly polarizable bodies it may be possible to neglect 
higher-order terms in the Born expansion and approxi- 
mate the scattering part of the Green tensor by the term 
linear in the susceptibility, and the decay rate accord- 
ingly. 

Using the real-cavity model of the local-field correc- 
tion, we have applied the theory to the problem of the 
local-field correction to the spontaneous-decay rate of an 
excited atom embedded in a dispersing and absorbing 
dielectric body of arbitrary size and shape. In this way, 
we have derived a rate formula, which, within the linear 
Born approximation, applies to atoms in arbitrary di- 
electric bodies, and have given explicit conditions of its 
validity. To illustrate the results, we have considered the 
case of an atom at an arbitrary position inside spherical 
body in more detail. 

It has surprisingly turned out that the rate formula 
found in linear Born approximation agrees, to linear or- 
der in the susceptibility, with a rate formula suggested 
by Tomas [l3| from his study of the real-cavity model 
in the special and analytically solvable case of a spheri- 
cally symmetric system. We have then shown that this 
quite general formula indeed applies to atoms in dielec- 
tric bodies of arbitrary sizes and shapes, provided that 
the atoms are not in the very vicinity of the surfaces of 
the bodies. So it can be shown that the scattering part 
of the Green tensor that enters the basic formula for the 
decay rate can always be decomposed into a term that 
only depends on the properties of the local environment 
of the guest atom and the local-field corrected scattering 
part of the Green tensor of the host body without the 
guest atom, where the correction simply appears in form 
of a factor. In particular, the former term can be sub- 
stantially determined by the absorptive properties of the 
host body, thereby giving rise to a shift of the decay rate. 

Finally, we note that in the same spirit as in the treat- 
ment of the spontaneous decay, the Born expansion can 
also be employed in studying other phenomena of the 
atom-field interaction in realistic systems whose Green 
tensors are not known or analytically too involved. Typ- 
ical examples may be the Casimir-Polder interaction of 
a ground-state atom with an inhomogeneous body |27| . 
the resonance fluorescence of an atom near such a body, 
or the resonant energy transfer between atoms embedded 
in realistic bodies. 
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